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This paper is devoted to a Radon-type transform arising in a version of Photoacoustic Tomography 
that uses integrating circular detectors. We show that the transform can be decomposed into the 
spherical Radon transform and the Funk-Minkowski transform. 

An inversion formula and a range description are obtained by combining existing inversion for¬ 
mulas and range descriptions for the spherical Radon transform and the Funk-Minkowski transform. 
Numerical simulations are performed to demonstrate our proposed algorithm. 
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Introduction 


Hybrid biomedical imaging modalities have been an active area of research lately, com¬ 
bining different physical signals in order to utilize their advantages to enhance images. 
Photoacoustic Tomography (PAT), which incorporates ultrasound and optical or radio¬ 
frequency electomagnetic waves, is one of the most successful examples of such a com¬ 
bination. While pure ultrasound imaging typically leads to high resolution images, the 
contrast between cancerous and healthy tissue is rather low. On the other hand, optical 
or radio-frequency electromagnetic imaging can provide a significant contrast; however, 
it has low resolution. The photoacoustic effect, discovered by A.G. Bell [^, allows one to 
take advantage of the strengths of the pure optical and ultrasound imaging, and it forms 
the basis of PAT. 

In PAT, an object of interest is irradiated by pulsed non-ionizing electromagnetic en¬ 
ergy. Due to the photoacoustic effect, an acoustic wave dependent on the electromagnetic 
absorption properties of the object is generated This wave is then measured by 

ultrasound detectors placed outside the object. The internal photoacoutic sources are 
reconstructed from the measurements to produce a 3-dimensional image (tomogram). 
Irradiated cancerous cells, in particular, are displayed with high contrast. This is due to 
the fact that such cells absorb several times more electromagnetic energy than healthy 
tissue, and electromagnetic deposition is proportional to the strength of the generated 
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acoustic wave. Since the initially generated acoustic wave contains diagnostic informa¬ 
tion, one of the classical mathematical problems of PAT is the recovery of this initial 
acoustic pressure held from ultrasound measurements made outside the object. 

The initial approach to detecting ultrasound signals in PAT has been to use small 
piezoelectric transducers. As these transducers mimic point-like measurements, recon¬ 
struction algorithms produce images with a spatial resolution limited by the size of the 
transducers. Another drawback of such detectors is the difficulty manufacturing small 
transducers with high sensitivity. To overcome these dehciencies, various other types of 
acoustic detectors have been introduced, e.g. linear, planar, cylindrical and circular de¬ 
tectors [i,[il,[il,[Il,[2i,[33. These detectors are modeled as measuring the integrals of 
the pressure over the shape of the detector. 

There are works discussing PAT with circular integrating detectors 124 
Some works 


M,\ 


35l-l37|]. 


37l | have considered PAT with circular integrating detectors on 


cylindrical geometry: a stack of parallel circular detectors whose centers lie on a cylinder. 
It was shown that the data from PAT with circular detectors is the solution of a certain 
initial value problem, and this fact was used to reduce the problem of recovery of the 
initial pressure to that of inversion of the circular Radon transform in The cases 

where the centers of the detector circles are located on a plane or a sphere have been 
discussed in [1^. In the present work we study PAT with circular integrating detectors 

small objects [35|. We note 
[35(. However, the geometry 


in a spherical geometry, which can be useful in investigating 
that a type of spherical geometry was also discussed in 
considered here is different, as only detectors with the same radius are required, while 
various sizes of circular detectors are needed in [35|]. Also, our approach is to define a 
Radon-type transform arising in this version of PAT, and we show that this transform can 
be decomposed into the spherical Radon transform and the Minkowski-Funk transform, 
both of which are well studied. 

This paper is organized as follows. Section [2] is devoted to the construction of a Radon- 
type transform arising in PAT with circular integrating detectors. In Section [3l we show 
that the Radon-type transform is the composition of the Minkowski-Funk transform and 
the spherical Radon transform with centers on a sphere, and we present the inversion 
formula using this fact. In Section 01 we describe the range of this Radon-type transform, 
and we give necessary and sufficient conditions for a function to arise as data measured 
by the circular detectors in the spherical geometry. Section [5] is concerned with numerical 
implementation. Section [6] contains some final remarks. The appendices provide proofs 
of some technical statements. 


2. Photoacoustic tomography with circular integrating detectors 

In PAT, the acoustic pressure p{x, t) satisfies the following initial value problem: 

d?p(x,t) = A^p(x,t) {x,t) G X (0, oo), 

p{x,0) = f{x) X = (xi,X2,3:3) G M^, (1) 

dtp{x, 0) = 0 X G 

Here we assume that the sound speed is equal to 1 everywhere, including the interior of 
the object. The goal of PAT is to recover the initial pressure / from measurements of p 
outside the support of /. 

Throughout this article, it is assumed that the initial pressure field / is smooth and 
supported in i?(0, r^et) - the ball in with radius centered at the origin. We also 
assume that the acoustic signals are measured by a set of circular detectors centered at 
the origin with radius - that is, the detector circles are great circles on the sphere 
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Figure 1. The great circle on the sphere dB{0, Vdet) of radius in the plane perpendicular to 9. 

dB{0,rdet) (see Figure [2|. 

The data P{6, t), for (0, t) £ S'^ x (0, oo) measured by the circular detector lying in the 
plane perpendicular to 9 (see Figure [2]) can be expressed as 

-P(0,t) = ^ / Pi'^deta,t)dS{a), 

§-a=0 


where dS is the measure on the unit circle S^. 

Everywhere in the text we will use the hat notation to indicate a unit vector in 2D or 
3D, and dS to denote the measure either on the circle or on the sphere. The dimension 
will be clear from the context. 

Using Kirchhoff’s formula for the solution of m (e-g- 0 ): 




p{x, t) = dt 


1 


4vrt J 

y dB{x,t) 


f0)dS{l3) , 


we can represent the measurements P{9,t) as follows: 

/ 


Pie,t) = 


1 


dt 


1 

dvrt 




§-a=0 \ dB(rdata,t) 


t 


fmdSW)j dS{a) 
f{rdetd + t^)dS0)dS{a) \ , 


\ e-a=Q 


( 2 ) 


where 5^ is the unit sphere in R^. Let us define the new Radon-type transform 

T^pf{d,t)= J j f{rdetd + t$)dS0)dS{a). 

e-&=o 
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Then, ([2]) reads as 


p{9,t) = ^dt{tnpf{o,t)). (3) 

If / is odd, i.e, f{x) = —f{—x), then TZpf is equal to zero. We thus assume that / is 
even. This entails no loss of generality, as any function compactly supported in the upper 
hemisphere of B{0,rdet) can be extended to an even function with support in B{0,rdet)- 
In what follows, TZp is seen to be the composition of the spherical Radon transform with 
centers on dB(0,rdet) and the Minkowski-Funk transform. This observation will allow 
us to recover / from the measurements P as well as to describe necessary and sufficient 
conditions for a given function to arise as data measured by the circular integrating 
detectors. 


3. Reconstruction 

The following definitions will be used in the rest of this paper. 

Definition I The Minkowski-Funk transform F maps a locally integrable function cf 
defined on S'^ x [0, oo) such that (j){a,t) = 4>{—a,t) into 

F(l){9,t)= j 4>{a,t)dS{a), /or (0, f) G 5^ x [0, oo). 

9-a=0 


Definition 2 The spherical Radon transform Rs maps a locally integrable function f 
on into 

Rsf{a, 0 = y firdeta + tp)dS0), for {a, f) G 5^ x [0, oo). 


Now it is easily seen that the following representation of TZpf holds. 

Theorem 3 For any f G C'°°(M^) compactly supported in B{0,rdet) with f{x) = f{—x), 
we have TZpf {9,t) = F{Rsf){9,t). 


Thus, in order to reconstruct the initial pressure /, we will employ known inversion 
formulas for the Minkowski-Funk and the spherical Radon transforms. A number of 
inversion formulas for these transforms have been derived, e.g. some of those for the 
spherical Radon transform can be found in lol. 1^. 291 and those for the Minkowski- 

Funk transform are studied in dJ, 1^1 181 . 1281] . We choose the inversion formula for the 
spherical Radon transform given in ajlSj and present the relation between the Minkowski- 
Funk transform and the regular Radon transform in the following theorem: 

Theorem 4 Let 4> be a continuous and even function on S'^, i.e., (/(0) = (l){—9). If 


Fm = 


(j){a)dS{a), 


e-a=o 


for 9 G 52 , 
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then, if 0 ^ ( 0 , 0 , 1 ) and 9 / ( 0 , 0 ,— 1 ), the following relation holds, 




where 9 = {9', 9^) = ( 01 , 02 ,^ 3 ) G S‘^,a = (a', 03 ) € <l>(a'/a 3 ) = 2(l){a)a\, and 


X2)](w, s) = / + z2t!;-*-)dz/, /or (a), s) G 5 ^ x M. 


The proof of this Theorem is provided in Appendix A. Our main result is given in the 
theorem below. 

Theorem 5 Let f G C'°°(M^) be an even function with compact support in B{Q,rdet)- 
Then, if G{uj, s,t) = (1 + s‘^)~ 2 'Rpf [u/y/l + s'^, —s/y/l + s'^, the following relation 
holds: 


fix) = 



dsdttdttG{u], s, t)\t=\r^,,&-x\ ds dS{u:) 


a^u ■ a' — sa\ 


\rdet6i - x\ 


dS{a) 


(4) 


Proof. From Theorem 01 we have 

F'(l)i9, t) = X2, t)] (^1^, “ ’ 


(5) 


where 9 = {9',9^), a = (a', 03) G 5^, <l>(a'/Q;3, t) = 2(j){di,t)a\ and R is the 2D Radon 
transform of with respect to the first two variables. That is. 


R[^{xi,X2,t)]{uj 


= y 


sch + vCj^,t)dv, 


for (w, s) G 5 ^ X R. 


After changing variables in ([5]) we obtain the relation 

1 


R[^{xi,X 2 ,t)]{Lb,S,t) = 


n/I + 




OJ 


—s 


Vl + ' Vl + 


t . 


( 6 ) 


In order to recover <hjWe use the filtered backprojection formula for the inversion of the 
2D Radon transform 


^>(xi,X2,t) =^iR*HdsR^)ixi,X2,t) 
47r 


1 


p.v. 


ds'iR<^>)iuj,s',t) 


s — s 


dSicb), 


(7) 


^For completeness, we verify the application of the filtered backprojection formula to the function <I> in Appendix 

B. 
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where R^u{xi,X 2 ) ■= / tt(a), (xi, X 2 ) •tD)hS'(w) is the backprojection operator and 
1 f u(s') 

Hu(s) := p.v.— / - - ds' is the Hilbert transform. 

Jr s-s' 

Applying ©to ([6]) leads to the inversion formula for the Minkowski-Funk transform: 


(p{a,t) =2 t) 


1 


R* I Hd, 


Sttoo V 




—s 


vi + Vvi + ' Vi + 


a 


as 


( 8 ) 


Now let = {Rsf){a,t). Theorem [3] implies that Ftp = TZpf. Hence, we obtain the 

following expression for the spherical Radon transform of /: 


Rsf{a,t) =^R#(Hd^ 


Svral \ 


1 


Vi +' 


K^p/) 




—s 


Vi + ’ Vi + 




R*Hds[G{u},s,t)] (—,t] 

ovrag \ as / 


a 

— 

as 


(9) 


In order to recover /, we apply the inversion formula for Rs given in 


/w = -^/ 

5= 


rdet f dttdttRsf{a,t) 


t=\rdet&-x\ 


dS{a). 


( 10 ) 


Combining two equations ([9]) and ([TO]) gives us 


fix) = I ^R* [Hds dttdttG{(b,s,t)] 


dSja) 

t=\rdet&-x\ \xdeta - x\ 


which is equivalent to ©. □ 

Corollary 6 Let f G be an even function with compact support in B{0,rdet)- 

Then, ifg{u,s,t) = {1 + s^)~ 2 P / y/\ -\- — s / \/l + s"^ ,t], the following relation holds 


fix) 



dsdttgiu), s, t)\t=\r^^,a-x\ ds dSiu) 


a^u} ■ a' — sal 


\rdeta - x\ 


dSia). 


Indeed, since Pi0,t) = (Svr^) ^dttlZpfi9,t) (see ©) we can write 

1 ^ _ dttGiui,s,t) 

{VTT^'VTT^'J Stt^VIT^’ 

where Giix, s, t) is the function dehned in the preceding theorem. Then the result follows 
immediately from Theorem [5j 
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4. Range description 


In this section, we describe the range of T^p and necessary and sufficient conditions for 
a function P to arise as data measured by the circular detectors. 

As mentioned before, the Minkowski-Funk transform and the spherical Radon trans¬ 
form are well studied, so their range descriptions have already been discussed (see [H- 
it was proved that the Minkowski-Funk transform is a (continuous) 
bijection from to itself, where = {(/>£ : 4>{9) = 

In 0, M. Agranovsky, P. Kuchment, and E. T. Quinto provided four different types 
of complete range descriptions for the spherical mean transform. In particular, one of 
the range descriptions states that necessary and sufficient conditions for the function 
g £ C^{S‘^ X [0,2rdet]) to be representable as g = Rgf for some / £ C“(R(0, r^et)), is 
that for any integers m,l, the (m,/)-th order spherical harmonic term (Hs_ig')m,z(A) of 

H|_i 5 (( 0 , A) vanishes at non-zeros of the Bessel function Jrn+n/ 2 -i(A), where 


H._i5(0,A) = I g{e,t)j^.i{Xt)p-^dt. 

0 

Here jnit) is the spherical Bessel function: 




where Jn(t) is the Bessel function of the first kind of order n and T(t) is the Gamma 
function. 

Combining these two range descriptions, we have the following range description for 

TZp: 

Theorem 7 Let g £ x [0,2rrfet]). The following conditions are necessary and 

sufficient for the function g to he representable as g = TZpf for some f £ C'“(R(0, r^et)).' 

• g is even in 9 £ S'^, i.e., g{9,t) = g{—9,t). 

• For any positive integer n, let 


OO 

H^„i5(0;A) = J g{9,t)j^.,iXt)e-^dt. 

0 

Then for any integers m,l, the {m,l)-th order spherical harmonic term 
(H|_ig()m,z(A) of H|i_i( 5 r( 0 , A) vanishes at non-zeros of the Bessel function 
Tm.+nj'l—X (A) ■ 

Proof. It is enough to show that the (m, l)-th order spherical harmonic term 
(H:|_ig')m,z(A) of H|_i 5 ( 6 *,A) vanishes at non-zeros of the Bessel function Jm+n/ 2 -iW 
if and only if the (m,/)-th order spherical harmonic term (Hi_i[E“^ 5 ])m,z(A) of 

H 2 _i[E“^g(]( 0 , A) vanishes at non-zeros of the Bessel function Jm+n/ 2 -i(A), where F~^ 
is the inversion of the Minkowski-Funk transform. 
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In (I^ . the inversion formula for the Funk-Minkowski transform is presented as follows: 


0 52 


u=l 


where <5 is the Dirac delta function. Therefore, F ^g{a, t) can be represented as 


1 


0 5= 


u=l 


Consider the (m, ^)-th order spherical harmonic term ^fi'])m,z(A): 


(H|_i[F-^5])m,z(A) = J U^.i[F-^g]{a,X)Yr{a)dSia) 

5= 

u 

1 

dvr^ 

1 

dvr^ 


52 


T / / A)i(« . 6 - 

0 52 


d 


Yr{a)dS{a) 

U=1 

dS{9)dv 


j Hn_i5(0,A) J 

k 0 52 52 


u=l 

( 11 ) 


where in the last line, we used the Fubini-Tonelli Theorem. 

To compute the inner integral with respect to a of (El), we need the Funk-Hecke 
theorem la, 


3l[ |: if \h{t)\dt < 00 , 


f h{e ■ u)Y{^{Qj)dS{uj) = c(n, m)Yi^{e), 

;;... 

-1 


Here > —1/2 is the Gegenbauer polynomial of degree m and Cm{t) = Pm{t) is 

the Legendre polynomial of degree m. Applying the Funk-Hecke theorem to the inner 
integral of ED, the (m,/)-th order spherical harmonic term (Hn-ilT ^fl'])m,z(A) is equal 
to 


1 


d 

du 


Hn_i5(0,A)yr(0)rfS(0) 


0 52 


{v? — 


dv 


U=1 
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Thus, we have 


(Hn_l[F-l5])rr^,KA) 


2tx 


d fPmiVl^) 
du J {v? - U2)V2 

l 0 


u=l 


( 12 ) 


Since for any integer m, there exists a positive number e such that PmW^ — v‘^){u^ — 
^ 2 ^-i /2 jg 2 ;ero on the interval [1 — e, 1], — v‘^){u^ — v‘^)~^/‘^dv is 

not constant on [1 — e, 1]. Therefore, from (fT^ . we have = 0 if and only 

if (Hn_l[F-l5])rr^,^(A)=0. ' ’ □ 


Since P{0,t) = {8TT‘^)~^dttTlpf{6,t), we have the following necessary and sufficient 
conditions for a function P to be data measured by the circular detector from the initial 
pressure field. 

Theorem 8 Let P G C^{S‘^ x [0,2r(iet])- The following conditions are necessary and 
sufficient for the funetion P to he data measured by the eircular detectors from the initial 
pressure field f G C^{B{0,rdet)): 

• P is even in 9 £ 3“^, i.e., P{9,t) = P{—9,t). 

• The integral of P with respeet to t from 0 to 2rdet is equal to zero. 

• For any positive integer n, let 


CO t 

H5_iP(0;a) = j J Pie,T)j^.,{xt)r-^drdt. 

0 0 

Then for any integers m,l, the {m,l)-th order spherieal harmonic term 
(Hi_iP)m,z(A) of H|._iP(0,A) vanishes at non-zeros of the Bessel function 

Jm+nl2—l (A) • 

Proof. Let us define a function g on 3“^ x [0, 2rrfez] by 

t 

90P) = P0,T)dT. 

0 

Then g is infinitely differentiable with respect to t away from zero, g{6, 0) = 0 is equivalent 
to P{9,0) = 0 by L’Hopital’s rule, and g{9,2rdet) = 0 is equivalent to 


2 ^" det 

J P{6, T)dT = 0. 

0 

It is sufficient to show that g is infinitely differentiable and continuous at 0. This will be 
accomplished in the following two lemmas, which are proved in Appendix A. 

Lemma 9 Let 4> £ C'°°([0, e]) for some e > 0 satisfy (j){0) = 0, and let n be a positive 
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integer. For any positive integer k <n, we have 

J ^{r)dt as t ^ 0. 

0 

Lemma [ 9 ] is used to prove: 

Lemma 10 Let (p £ C'°°([0,e]) satisfy 0(0) = 0, and let h{t) = 4>{T)dT/t. Then 

h £ C°°([0,e]) and h{0) = 0. 

Theorem [7] now immediately follows by applying Lemma [10] to g{9,t) considered as a 
function of t for fixed 9. □ 


5. Numerical Experiments 

In this section, we show the results of numerical experiments reconstructing the initial 
pressure field / from the data P{9,t) = -^dt (tTZpf{9,t)^. We set rdet = 1 and assume 

that P{9,t) is given on a spherical grid in 9 on the unit sphere and a uniform grid in t 
on the interval [0,2]. 

We follow a method of reconstruction using the formula in Theorem O 

(1) Compute t) for each fixed t to get the values of p{9, t) on S‘^ x [0, 2] 

(2) Invert the wave operator to get f{x) 

Our phantoms are taken to be either the sum of multiples of characteristic functions 
of balls 


Xsh(^) 


c if |x — xoj < r 
0 if |x — xoj > r 


or the sum of C^-functions of the form 


(13) 


Xsm(x) = op (|T- aToj^) , 


(14) 


where p is defined by 


T (s^) = 


1-3 


s — r 
R — r 


+ 2 


s — r 
R — r 


if jsj < r 
if r < jsj < i? 
if |s| > R 


and 0 < r < R. Our algorithm reconstructs the even part of the phantom, so we take 
phantoms supported in the upper half of the unit ball where X 3 > 0 . 

The acoustic pressure generated by Xsh is 


p{x,t) 



if \x — xo\ — t < f 
ifjx — xoj— 
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and the acoustic pressure generated by Xsm is given by 


(see [13, Theorem 5.2]). In order to get the data P{9,t) measured by our detectors, we 
use a quadrature on each detector circle at each grid point in t, with the trapezoid rule 
and 100 points. 



(a) Phantom (b) Cross-section of the reconstruction in the X 2 Xs- 

plane 




(c) Slice of the reconstruction and phantom along (d) Slice of the reconstruction and phantom along 
the a: 3 -axis the line xi =0, xs = 0.5 


Figure 2. Reconstructions of a sharp phantom with centers (0,—0.6,0.4), (0,0,0.5), (0,0.55,0.5), and 
(0,0.55,0.15), radii 0.2, 0.3, 0.15, and 0.1 and values 0.6, 1, 1.1, and 0.8, respectively 


The Funk-Minkowski transform is inverted at each time step by using the inverse two- 
dimensional Radon transform as prescribed by Theorem [4l A direct inversion of F using 
Theorem m would involve dividing by which could cause instability near the aia 2 - 
plane. To avoid this, one possible remedy is to cyclically permute 9i, 62 , and 6*3 in ([5]) 
and sum the results to obtain 

2p((a, t)((a]^ cx2 Trr3) — 2p((a, t). 

In fact, we can further reduce the contribution from points near the a* = 0 plane (where 
ai is the coordinate in the denominator of ([5]) after a cyclic permutation of 9i, @ 2 , and 
@ 3 ). To achieve this we multiply the result of ([5]) by and sum the results: 

2 p{a, t) (ai+^ -h 02 +^ + ^ 3 +^). 

The function p{a,t) is then easily obtained by dividing by 2 ( 0 ^^^ -|- In 

our reconstructions, we use A: = 2 in this formula. To get the initial pressure field from 
p{a,t), we use (fT 0 |) . 

In our experiments, the function p{9, t) is reconstructed on a spherical grid of 50 polar 
angles and 200 azimuthal angles, for a grid of 50 values of t on the interval [0,2]. The 
polar angles range from [7r/25,7r/2]. We start our polar grid slightly away from 0 because 
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0.8 



(a) Phantom 




0.8 



(b) Cross-section of the reconstruction in the X 2 X^- (c) Cross-section of the reconstruction in the 3 : 23 : 3 - 
plane from exact data plane from data with 20 % uniform random noise 




(d) Slice of the reconstruction and phantom along (e) Slice of the reconstruction and phantom along 
the 3 : 3 -axis from exact data the 3 : 3 -axis from data with 20 % uniform random 

noise 




(f) Slice of the reconstruction and phantom along (g) Slice of the reconstruction and phantom along 
the line xi =0, 3:3 = 0.5 from exact data the line xi = 0, 3:3 = 0.5 from data with 20% 

uniform random noise 


Figure 3. Reconstructions of a smooth phantom with centers (0,—0.6,0.4), (0,0,0.5), (0,0.55,0.5), and 
(0,0.55,0.15), outer radii 0.2, 0.3, 0.15, and 0.1, inner radii 0.1, 0.15, 0.075, and 0.05, and central values 0.6, 
1, 1.1, and 0.8, respectively 


when we invert the Funk-Minkowski transform using Theorem U small values of the 
polar angle correspond to s values near infinity, and having to account for large values of 
s makes implementation of the inverse Radon transform difficult. It is also not necessary 
to take polar angles past 7r/2 because we are reconstructing the even part of /. Our final 
reconstructions are calculated on an 80 x 80 x 80 grid in x. 

Figure [2] shows the reconstruction of a phantom comprising four sharp balls of the form 
m from exact data. A reconstruction of a smooth phantom comprising four functions 
of the form (jl4l) from exact and noisy data is shown in Figure [3l 
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6. Conclusion 


We studied a Radon-type transform arising in PAT with circular detectors and showed 
that it can be decomposed into the Minkowski-Funk transform and the spherical Radon 
transform. An inversion formula and range conditions were provided. In addition, numer¬ 
ical experiments demonstrated the feasibility of the reconstruction formula and showed 
good reconstructions from noisy data. 
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Appendix A. Proofs of Theorem |4] and Lemmas [9] and 1101 

Theorem m Let 4> be a continuous and even function on , i.e., 4>{6) = If 

F<p{e)= J ^{a)dS{a), forOeS^, 

9-a=0 


then, if 6 ( 0 , 0 , 1 ) and 6 / ( 0 , 0 , — 1 ), the following relation holds, 

1 / 0' 9^ 


where 9 = {9', 9^) = ( 0 i, 6*2,6*3) G S‘^,a = (o', 03) G S"^, ^(a'/as) = 2(l){a)a^, and 


(Al) 


-R[<i>(xi, X2)](w, s) = / ^{sLJ + u 6 j'^)diy, /or (a), s) G x M. 


Proof. By definition, F(p can be written as 


2 tt 


Fm = / 4 >{u cos t + V sin t)dt, 


where 


U = {UI,U2,U5) 


-^(- 02 , 6 * 1 ,0) G S'^ and v = {vi,vi,V 3 ) 

P I 


1^(0103,0203,-10?) G 


Note that u,v and 0 are orthonormal, so that ucost + hsint, for t G [0,27r] is a 
parametrization of the great circle orthogonal to 0 . 

By the evenness of cf and the dehnition of <h, we have 



K 


where in the last line, we changed the variables cott —)• t. Here u' = (— 02,0i)/|0^| and 
v' = 0'/?|0'p + |0?/0|. Thus Ff{9) is the integral of <1> over the line with direction 
u'/v-i and passing through v' jvj,. We can represent Fcf as 


F(j){9) = / $ { + 


v' luO 


\V'\ V3 


Ug ‘^dt = 


^ i u't + 


v' \v'\ 


|v| 


|u3| ^dt 


= |u3|-'i?^>{^,^), 

' 1'^ I Vs 
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where in the second line, we changed the variables t/v^, t. The following identities 
complete our proof; 


v' 9' |r;'l 03 , I ,_i 1 

W\^wy ~wy ^ Wy 


□ 

Lemma [9] Let cj) e C'°°([0,e]) for some e > 0 satisfy 4>{0) = 0, and let n be a positive 
integer. For any positive integer k <n, we have 


Z 

J 

0 


(t-r) 


n+l 


t 


dr —)• 0 as t —)• 0. 


Proof. It is enough to show that for any positive integer m, 

t 


j 


{t - r)™' 


tm 


dr ^ 0 as t —^ 0, for any m' >m 


(A2) 


0 


since 


dy 





where Ck,j are some constants depending on k and j. To show (IA2I) . it is enough that for 
any positive integer m, 



tm 


-dr —)■ 0 


0 


as t —y 0, 


for any m' >m 


(A3) 


since 


{t-rr' = Y, 

j=0 


where m'Cj is the number of ^-combination from a given set of m' elements. Applying 
L’Hopital’s rule to ()A3D . we have 






as t —>■ 0. 


□ 

Lemma [To] Let 4> e C'°°([0,e]) satisfy ^(0) = 0, and let h{t) = 4>{T)dT/t. Then 

h e C'°°([0, ej) and h{0) = 0. 

Proof. The smoothness of h is obvious except at t = 0. In order to prove smoothness 
at t = 0, we will use mathematical induction. First of all, /i(0) is equal to (/)(0), which 
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is equal to zero, by L’Hopital’s rule. Suppose ^^(0) = (j)^^ ^)(0)/n and ^\t) is 
continuous at 0. Consider the Taylor expansion of (j){T)dT at t = 0. Then for any re, 
we have 


j (j){T)dT =t(j){0) + ^(/)'(0) + H- 

0 


fU 

+ 0(-i)(O) + 

re! 


fin+l) 

(re + 1)! 



^(n+l)(T) 

(re + 1)! ^ 


— r)"'+^dr, 


where is the re-th derivative of (p evaluated at t = 0. Hence, we have 

t 

h{t) J (j){T)dT 
0 

=0(0) + |f0'(O) + ... + i;^0<'»(O) + 1 « - r)»«0r, 

0 

For any k < n, the /c-derivative of the last term (or the reminder term) is 


j 0^l)(r) (( - I 0<"+‘>(r) / 

(n + 1)! t ^ J {n + iy.\ t J 

0 0 


(A4) 



0 


since 5|''((t —= (t — r)"'”''^ ^"[(re + l)!/(re —-1-(—l)^'A:'!(t — 

for k' < n. Hence, for t near 0, we have 




l<^(n-i)(o) + ^<^W(o)+ f 

re re + 1 ,/ (?^ + Ij! 


(+^)" 


and 






dr. 


By Lemma El we have 


hH (t) ^ (0) as t ^ 0. 

re + 1 
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Hence we have 


/i(")(0) = lim- 




n + 1 



(n + 1)! * 


(i-r) 


n+1 


dr 


n + 1 ’ 


where in the second line, we used L’Hopital’s rule, (IA4p . and Lemma[9l This means that 
h is n-th differentiable at 0 and is continuous at 0. □ 


Appendix B. Verification of ([7]) 


It is natural to apply an inversion formula for the regular Radon transform to obtain 
the inversion formula for the Minkowski-Funk transform when the relation (lAip is given. 
However, the filtered backprojection formulcd is typically stated for functions in the 
Schwartz space (see Em, [ 2 ^). For completeness, we verify here that the filtered back- 
projection formula can still be applied to Q. 

From now on, we omit t, an unnecessary variable here. From ^{a'/u'^) = 2(j){a)a'l for 
d = (a', as) G 5^, we notice that 


$(xi,X2) = (p 


(X1,X2, 1) 
y/l + x'^ + X2 


1 

1 + xl + X 2 


is a smooth and L^-function, i.e.. 


/ \^{xi,X2)\‘^dxidx2 < 00 . 

R2 


By the projection slice theorem, we have 4>((Tch) = R^{uj,a) for a Schwartz function d>, 
where ^ is the 2D Fourier transform of and is the ID Fourier transform of R^. By 
an approximation argument, the projection slice theorem still holds for an L^-function 
<k. Since the Fourier transform is a unitary isomorphism on the L^-space, the set of 


^In [ifll Thoerem 3.1 in Chapter 1], the function <I> can be recovered from the Radon transform by the following 
inversion formula 

^{xi,X2) = ^ [ \{xi,X2) - iyi,y2}\~^R*R^iyi,y2)dyidy2 
47r J 

provided $(a:i,X 2 ) = 0(|(a:i,ai2)|~‘^) for some > 1. However, the filtered backprojection formula we want to 
apply is $( 011 , 0 : 2 ) = ( 477 )~^HR^. 
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L^-functions, we have 


^{xi,X 2 ) 


[Z7r)^ R-^oo J 

\{ii,i2)\<R 

. lim [ [ 

{ZnyR^aoJ J V , I 

51 a<R 

—^ lim [ [ isgn(a)dJ^(uj,a)e^'^^^^’^^'^'^dad6j, 
[2ttYr-^coJ j 

51 a<R 


(Bl) 

where in the second line, we changed the variables (^1,^2) {<x,dj) £ [0,00) x S^. Since 

F(j) is smooth and bounded, i?<h(a),s) is a smooth and L^-function and so is dgR^- We 
notice that Hilbert transform Hu is also a unitary isomorphism on the L^-space (by the 
Titchmarsh theorem (^ . Chapter 5]). Therefore, (IBip becomes (j7|). 
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